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PREFACE 


Meteorological  studies  suggest  that  technologically  feasible  oper¬ 
ations  might  trigger  substantial  changes  in  the  climate  over  broad  re¬ 
gions  of  the  globe.  Depending  on  their  character,  location,  and  scale, 
these  changes  might  be  both  deleterious  and  irreversible.  If  a  foreign 
power  were  to  bring  about  such  perturbations  either  overtly  or  covertly, 
either  maliciously  or  heedlessly,  the  results  might  be  seriously  detri¬ 
mental  to  the  security  and  welfare  of  this  country.  So  that  the  United 
States  may  react  rationally  and  effectively  to  any  such  actions,  it  is 
essential  that  we  have  the  capability  to:  (1)  evaluate  all  consequen¬ 
ces  of  a  variety  of  possible  actions  that  might  modify  the  dim  e, 

(2)  detect  trends  in  the  plobal  circulation  that  presage  changes  in  the 
climate,  either  "atural  or  artificial,  and  (3)  determine,  if  possible, 
upeann  to  counter  potentially  deleterious  climatic  changes.  Our  posses¬ 
sion  of  this  capability  would  make  incautious  experimentation  unneces¬ 
sary,  and  would  tend  to  deter  malicious  manipulation.  To  this  end,  the 
Advanced  Research  Projects  Agency  initiated  a  study  of  the  dynamics  of 
climate  to  evaluate  the  effect  on  climate  of  environmental  perturba¬ 
tions.  The  presont  Report  is  a  technical  contribution  to  this  larger 
study. 

The  ability  to  discern  changes  in  the  general  circulation  as  sim¬ 
ulated  by  a  numerical  model  has  been  questioned  because  of  the  inhurent 
unpredictability  of  the  atmosphere.  This  study  was  undertaken  to  dem¬ 
onstrate  that,  despite  the  errors  that  inevitably  occur  in  day-to-day 
predictions,  the  mean  values  are  sensitive  to  major  changes  in  the 
boundary  conditions.  A  set  of  calculations  has  been  made  using  the 
Hinta-Arakawa  2-level  model  to  determine  whether  changes  in  the  mean 
state  can  be  found  in  the  presence  of  error. 

Reports  on  related  work  are:  The  fien*  *  t  f  the  Arctic  tfmtn 

tf-ai  !tfi  bv  0.  Fletcher,  R-444-PR;  A  .'v-f’yv’wta- 

tim  of  t he  tVr-.Viv*  Atnr-erherio  ra! 

Model,  bv  W.  I..  dates  el  al.,  R-877-ARPA;  and  iViii  Tuti'idi tu  .’’to dies, 
I:  VolanKie  ./u.u  —  A  >*t :!„•■»!  .‘ur:vy ,  bv  l),  D 1  emend )  i  an, 

R-886-ARPA, 
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SUMMARY 

The  growth  of  small  errors  in  numerical  models  of  the  atmospheric 
circulation  destroys  the  detailed  predictive  capability  of  those  models 
within  a  few  days.  Despite  the  failure  of  the  models  to  produce  accu¬ 
rate  local  predictions,  it  was  hypothesized  that  a  change  in  the  equator- 
to-pole  temperature  gradient  would  produce  discemable  effects  in  average 
conditions.  This  Keport  presents  the  results  of  an  experiment  to  test 
this  hypothesis. 

The  Mintz-Arakawa  model  was  started  with  a  standard  set  of  initial 
conditions.’  and  was  run  60  days.  The  experiment  was  then  replicated 
twice,  with  two  independent  sets  of  random  temperature  variations  super¬ 
posed  on  the  temperatures  at  the  •)  »  0.25  and  j  ■=  0,75  levels.  At 
this  point,  the  ice  of  the  Arctic  Ocean  was  replaced  with  water  at 
the  freezing  temperature.  Again  the  model  was  run,  starting  once 
with  the  standard  initial  conditions  (other  than  the  ice  removal) 
and  once  each  with  the  two  sets  of  temperature  "errors"  added. 

For  the  four  replications  with  temperature  "errors,'  the  detailed 
predictive  capability  was  lost  after  about  14  days.  Yet  an  analysis 
of  variance  applied  to  the  lost  30-day  zonal-average  values  of  the  three 
•eta  of  ice-in/ieo-out  runs  showed  many  significant  changes  in  the  gen¬ 
eral  circulation.  These  significant  changes  may  be  summarized  as  fol¬ 
lows: 

(1)  The  vertical  stability  over  the  Arctic  Basin  decreased 
markedly.  An  expected  result  of  this  decreased  stability 
(and  the  newly  available  surface  moisture)  would  he  an 
Increase  in  ptecipitation;  analysis  of  this  possibility 
will  he  deferred. 

(2)  The  excfSo  heat  in  the  Arctic  caused  a  slight  efflux  of 
mass  from  tk,e  area,  vhtcn  caused  tower  pressure  ever  the 
Arctic  basin. 

(3)  The  oceanic  heat  source  present  upon  removal  of  the  ice 
weakened  the  cold  core  vortex  over  the  Arctic  Basin. 
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(4)  The  geopotential  heights  at  high  levels  in  the  sub-Arctic 
were  increased  by  this  change  in  the  distribution  of  mass. 

(5)  The  midlatitude  westerlies  were  weakened  by  the  redistri¬ 
bution  of  mass. 

These  results  show  that  changes  in  boundary  conditions  in  the  model 
produce  detectable  changes  in  the  general  circulation  despite  the 
breakdown  of  predictability. 
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I.  INTRODUCTION 


In  a  recent  paper,  Lorenz  (1969)  raised  a  serious  question  regard¬ 
ing  the  ability  to  predict  the  future  state  of  fluid  systems.  He  states 

...that  certain  formally  deterministic  fluid 
systems  which  possess  many  scales  of  motion  are 
observational ly  indistinguishable  from  indeter¬ 
ministic  systems;  specifically,  that  two  states 
of  the  system  differing  initially  by  a  small  'ob¬ 
servational  error'  will  evolve  into  two  states 
differing  as  greatly  as  randomly  chosen  states 
of  the  system  within  a  finite  interval,  which 
cannot  be  lengthened  by  reducing  the  amplitude 
of  the  initial  error. 

Lorenz  go?®  on  to  say: 

It  is  found  that  each  scale  of  motion  possesses 
an  intrinsic  finite  range  of  predictability,  pro¬ 
vided  that  the  total  energy  of  the  system  does  not 
fall  off  tot*  rapidly  with  decreasing  wavelength . 

With  the  those n  values  of  the  constants,  'cumulus- 
seale'  Motions  cart  be  predicted  for  about  an  ,«our, 
'svnoptic-scale*  motion#  a  few  days,  and  the  larg¬ 
est  scale*  a  tew  weeks  In  advance. 

For  those  who  propose  using  siobal -ci  rets  la?  ion  models  to  investi¬ 
gate  purposeful  and  inadvertent  eiifcate.  fsodit  i  cat  ion  the  mol  scat  tuns 
of  l.awo’i  thrMt  y»s  1  -*  s-  ser f indeed.  5*  or  tisaet  f  >■•«* .  .error# 
id  initial  condition*  would  tend  to  taa  v,  tespirirlv  the  desired  effect 
in  a  relatively  *m*» t  lit**,  there  are.  however,  Soar  ittp« riant  differ¬ 
ence;,  between  the  primitive  equation  eedet1*  and  the  idealised  svstea 
from  which  Lorens  readied  hi*  conclusion*.  These  occur  erirvifilv  in 
th*  form  of  beating  Wfitss  functions  and  dissipative  tern*,  f olh  of 
which  are  ttissing  in  the  lafenx  ferwilation. 

Ouf  Purpose  1*  t’o  in-'vst  irate  the  ee**iM  Htv  of  dele,  ting  rhaftees 
in  the  mean  condition®  du«  to  beeodarv  value  rba«er«  so«eitv  the  fact 
that  predict. •rh 5 1 .3 tv ,  «n  th*,  Loren*  sense,  ha*  disappeared  before  the 
period  of  averaging  j>-  cenplcle. 

In  the  absence  of  an  onaivsis  ei  error  growth  in  the  ps  eitive 
equations  of  motion.  We  loni  ior  at’  et  neaps  to  assess  the  IMeUtv 
of  current  global  ct  reolatirsi  taodels.  ihis  report  present*  the  results 


of  a  series  of  six  60 -day  integrations  of  the  Mintz-Arakawa  two-level 
global-circulation  model,  designed  to  test  the  hypothesis  that  the 
presence  of  sizable  changes  the  boundary  conditions  could  be  deter¬ 
mined  in  the  solutions.  Ihe  variables  considered  are  values  (zonally 
averaged  and  time  averaged)  from  the  last  30  days  of  each;  60-day  run, 
and  the  significance  of  the  changes  is  tested  by  classical  analysis- 
of-variance  technique.  Ihe  average  differences  between  two  arctic  ice 
conditions  are  then  discussed  in  terms  of  qualitative  hypotheses  that 
have  been  advanced  to  explain  the  circulation  changes. 
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II.  THE  EXPERIMENT 

It  should  be  understood  that  the  experiment  described  herein  is 
one  that  estimates  only  the  model 's  ability  to  separate  signal  (pur¬ 
poseful  changes  in  boundary  conditions)  from  noise  (random  changes  in 
initial  conditions).  No  claim  can  be  made  that  the  real  atmosphere 
would  behave  in  a  similar  way;  one  hopes  that  it  would,  but  this  is 
a  separate  problem  and  is  noc  addressed  here. 

The  experiment  is  easily  described  in  a  step-by--step  manner.  Our 
initial  conditions  were  taken  from  day  370  of  die  Mincz-Arakava  model 
as  run  by  UC1A.  From  this  standard  set  of  initial  and  boundary  condi¬ 
tions  for  die  model  (Gates  et  ai.,  1971),  an  initial  "control  run"  of 
sixty  days  was  made.  The  first  th’rty  davs  were  to  allow  transients 
to  settle  and  boundary  conditions  to  make  theii  effect  felt.  The  last 
thirty  days  was  the  period  over  wi  £  eh  the  results  were  averaged.  Next, 
this  run  was  repeated  exactly,  hut  with  a  small  random  et  rot  added  to 
the  initial  free-air  temperature  conditions  at  each  grid  point.  The 
distribution  from  which  this  random  error  was  drawn  had  zero  mean  and 
a  variance  of  1*C  (less  chon  the  variance  from  actual  measurements). 
These  additive  errors  say  be  viewed  as  reflecting  our  uncertainty 
about  global  initial  conditions  if  we  were  to  atu-espt  real  prediction 
with  the  model. 

After  completing  the  second  run-,  a  third  run  exactly  like  the  sec¬ 
ond  was  evade,  maintaining  the  sane  hound  arc  condition* .  but  with  a  cots* 
pletely  aew  set  of  random  noise  samples  drawn  ft or-  the  same  distribution 
a*  before.  W*  new  had  three  samples  with  the  »**»  boundary  conditions 
and  a-nost  the  sar*?  initial  conditions,  whether  tie  selut  ions'  as  a  func 
tics  of  tint1  would  hf  -suite  to  each  ether  or  « on  *  f  I  cant  1\  differ¬ 

ent  Was  unknown .  1W  *»mer*  |S  uM  els  the«c  «  '7  I  .«s\*  ire  <ens  Sdervd 
"close"  depended  >*.  tl;<.  •>».•*. c **s*r  ei  t;.«-  =.«-  *  ••  >  *  *»..s,le  described 

f  e  1  e»i . 

At  the  beginning  of  3*un  « ,  w  introduced  a  -oar.*-.-  it  the  outface 
boundarv  conditions,  discussed  below.  The  Mrs!  run  with  t?a ».*  new  bound 
arv  condition  was  -strain  cal  Sod  ti  centre’  fun  and  ,*d  no  additive  randov 
noise  cent  aeir.al 5nr  ti  r  initial  conditions,  ^uns  >  aid  b  do  haw-  > ho 


additive  noise  and  this  noise  is  identical  to  that  which  perturbed  the 
second  and  third  runs  with  the  original  boundary  condition. 
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III.  CHANGING  THE  INITIAL  CONDITIONS 


We  elected  to  perturb  the  initial  temperature  at  the  two  Interior 
levels  of  the  model  as  the  means  of  introducing  noise  into  the  initial 
conditions;  these  levels  are  identified  in  the  model  as  a  =  0.2S  and 
a  -  G.75,  and  are  approximately  the  400  and  S00  mb  surfaces,  respectively. 
The  temperature  values  at  each  grid  point  were  modified  by  a  number 
drawn  from  a  normal  distribution  having  zero  mean  and  a  variance  corre¬ 
sponding  at  most  to  the  variance  encountered  in  actual  temperature 
measurements.  Drawing  from  a  zero-mean  population  virtually  assured 
that  the  post-perturbation  global  mean  temperature  would  remain  un¬ 
changed  (given  sufficient  grid  points). 

The  choice  of  magnitude  of  the  variance  of  the  noise  was  more  dif¬ 
ficult.  If  a  large  variance  were  chosen,  the  solutions  might  diverge 
drastically,  leaving  us  to  defend  the  challenge  that  the  noise  could 
realistically  have  been  smaller  with  possibly  a  different  outcome.  If 
a  small  variance  were  chosen  we  would  have  been  subject  to  challenge 
for  not  making  a  significant  change.  Of  course,  if  very  small  pertur¬ 
bations  had  caused  very  different  solutions,  an  important  characteris¬ 
tic  of  the  model  would  have  been  exposed.  We  chose  the  variance  of  the 
temperature  perturbation  to  be  1°C.  This  is  in  accord  with  current  GARP 
data  specifications  and  has  been  used  in  other  predictability  experi¬ 
ments  (Jastrow  and  Haletn,  19  70,  Chamey,  1966). 

There  are  initial  errors  other  than  temperature  that  could  affect 
the  nature  of  the  solutions,  but  our  main  objective  was  to  ensure  only 
that  solutions  with  perturbations  were  different  from  control  runs  — 
not  to  study  the  nature  of  the  errors  as  such  .  In  a  later  section  we 
attempt  to  show  that  the  temperature  errors  introduced  did,  indeed, 
produce  the  type  of  unpredictability  discussed  by  Lorenz. 
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IV.  CHANGING  THE  BOUNDARY  CONDITIONS 


The  philosophy  involved  in  choosing  a  change  in  boundary  condition 
is  considerably  more  subtle  than  that  of  the  initial-condition  pertur¬ 
bation.  There  is  little  scientific  basis  or  prior  experience  to  draw 
on.  The  purpose  of  the  change  is  to  induce  a  substantial  change  in  the 
climate.  In  our  case,  ’’substantial"  could  mean  Just  large  enough  to 
be  perceived  through  the  noisy  initial  conditions,  or  it  could  mean 
large  enough  to  alter  the  human  environment  significantly  somewhere  on 
Earth,  possibly  distant  from  the  immediate  locale  of  the  change. 

Examples  of  major  changes  in  boundary  conditions  could  be  the  re¬ 
moval  of  all  or  a  portion  of  the  arctic  ice,  flooding  the  Sahara  Desert, 
or  radically  changing  the  evaporation  of,  say,  10  6  sq  ktn  of  low-latitude 
ocean.  There  was  no  experience  which  led  us  to  believe  with  certainty 
that  the  model  would  respond  to  any  of  these  changes  with  a  substantial 
climatological  change.  They  are  "major"  only  because  they  represent 
»ery  difficult  or  nearly  impossible  engineering  feats  for  man. 

There  was  no  compelling  reason  to  subject  the  system  to  such  a  se¬ 
vere  change.  An  equally  interesting  approach  would  have  been  to  try  as 
small  a  change  as  one  would  guess  was  capable  of  producing  an  observable 
effect.  Here  the  risk  was  that  the  effect  would  not  be  discernible. 

In  the  experiments  reported  here  we  chose  to  alter  the  surface¬ 
boundary  condition  by  removing  all  the  sea  ice  in  the  northern  hemi¬ 
sphere  and  replacing  it  with  water  at  *-l®C.  In  the  Mintz-Arakawa  model,* 
the  ocean's  surface  teuperature  is  prescribed,  so  once  set  at  -1°C  the 
arctic  water  temperature  remains  at  that  temperature. 


A  cooiplete  documentation  of  the  Mintz-Arakawa  model  appears  in 
Gates  et  al. ,  1971. 
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V.  PREDICTABILITY 


A  measure  often  used  to  estimate  predictability  is  the  root -mean- 
square  difference  between  a  control  run  and  its  perturbed  counterpart. 
Let  T^o^(t)  be  the  temperature  at  grid  point  (i,j),  sigma  level 
a  »  0.25,  0.75,  run  number  £  and  time  t.  The  run  number,  £  =  1,2,. ..,6, 
indexes  the  simulation  as  defined  in  Table  1. 


Table  1 

£- INDEXING  OF  EXPERIMENTS 

A  Experiment  Description 

1  Ice-in  control 

2  Ice-out  control 

3  Ice-in  perturbation  no.  1 

4  Ice-out  perturbation  no.  1 

5  Ice-in  perturbation  no.  2 

6  Ice-out  perturbation  no.  2 


Then  predictability  in  terms  of  temperature  is  measured  by 


ATo(t) 


(t)  -  T 


ijoB 


(t)  ] 


2|l/2 

f 


(1) 


where  B  «  3  or  5  if  a  «  1  and  B  ■  4  or  6  if  a  «  2.  Figures  1  to  4  show 

AT  (t)  for  o  e  0.25  and  0.75  and  a  =  1  and  2. 

0 

In  all  cases  the  value  of  AT^Ct)  dropped  to  about  0.5#C  in  two  days 
after  the  perturbations  were  introduced  (at  day  370).  This  is  most 
likely  due  to  the  smoothing  out  of  the  initial  random  errors  which  were 
introduced.  From  day  372  the  errors  increase  in  a  rather  regular  fash¬ 
ion  until  about  day  390.  From  then  until  the  end  of  the  run,  AT^Ct) 
oscillates  erratically  with  an  amplitude  of  about  2°C.  We  take  this  to 
indicate  that  the  systems  reached  states  where  the  predictability  has 
essentially  been  lost.  We  therefore  assume  that  any  statistically  sig¬ 
nificant  differences  between  means  over  the  Inst  thirty  davs  of  the 
runs  must  be  the  effect  of  changing  the  boundary  conditions. 


Day 


Fig.  1  —  Global  root-mean-square  temperature  deviation  as  a 
function  of  time:  ice  in;  o  “  0.75. 


--  Globa  1  root-mean-square  temperature  deviation 
function  of  time:  ice  in;  <?  •  0.25. 
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Day 


F 1  ^ .  3  —  Global  root- scan-squart*  temper  .itufv  leviaU-'o  .is  « 
function  of  t. tew :  ice  out;  ?  • 
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VI.  STATISTICAL  TESTS 

The  question  we  ask  is  whether  or  not  the  change  in  boundary  con¬ 
dition  produces  a  significant  change  in  some  climatic  variable  of  in¬ 
terest.  This  could  have  f ai led  to  occur  for  several  reasons:  (1)  the 
random  perturbations  in  initial  conditions  produced  inordinately  large 
changes  in  the  final  climatic  variables,  thus  totally  obscuring  the 
actual  effect  of  different  boundary  conditions,  (2)  the  boundary  con¬ 
dition  change  was  not  large  enough  to  produce  an  effect  in  the  allotted 
simulation  time;  or  (3)  the  change  in  boundary  condition  was  not  physi¬ 
cally  significant,  i.e.,  it  was  either  a  wrong  choice  for  affecting  fu¬ 
ture  climate,  or  the  model  did  not  properly  represent  the  physics.  A 
straightforward  analysis-of-variance  procedure  may  be  used  to  test  the 
hypotheses  that  the  boundary-condition  change  had  no  effect,  or  that 
the  additive  noise  had  no  effect. 

For  all  variables  which  are  averaged  over  the  last  30  days  and  over 
all  values  of  longitude  at  a  fixed  latitude  (the  "JO-dav  zonal  01030" 
the  statistical  procedure  reduces  to  a  simple  univariate  analysis  o i 
variance.  By  such  an  averaging  procedure  all  data  were  reduced,  in  our 
experiment,  to  a  2  *  1  array  of  scalar  variables  shown  in  Table  2. 


Table  2 

ARRAY  OF  SCALAR  VAR  1  ABU'S 


Unperturbed 

Perturbation  No.  1 

ehheehi 

Ice  In 

Y11 

Y12 

V13 

V 

Ice  Out 

Y21 

'22 

V2J 

■a 

YU 

*.> 

V*  3 

where  is  the  zonal ly  averaged  30-dav  mean  of  the  variable  of  interest. 

By  assuming  that  the  Y^  are  Independent,  normally  distributed  ran¬ 
dom  variables,  we  can  test  the  hypothesis  that  the  row  means  (Y  )  or  the 
column  means  (V  )  in  Table  2  come  from  the  same  distribution.  He. erring 
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to  Table  1  we  see  that  rejecting  the  null  hypothesis  that  all  row  means 
come  from  the  same  distribution  is  the  same  as  rejecting  the  hypothesis 
that  removing  the  arctic  ice  had  no  effect  on  the  physical  variable  un¬ 
der  observation.  Similarly,  rejecting  the  null  hypothesis  that  all 
column  means  come  from  the  same  distribution  is  the  same  as  rejecting 
the  hypothesis  that  perturbing  the  initial  conditions  had  no  effect  on 
the  variable. 

The  methods  available  allow  us  to  guard  only  against  type  I  errors. 
That  is,  when  the  null  hypothesis  is  true,  the  probability  of  rejecting 
it  can  be  made  as  small  as  we  wish. 

We  cannot,  however,  make  any  statement  about  the  probability  of  ac¬ 
cepting  the  null  hypothesis  when,  in  fact,  it  is  false  (a  Type  II  error). 
This  is  the  less  objectionable  of  the  two  types  of  error,  and  while  an 
estimate  of  its  probability  would  be  useful,  it  is  by  no  means  crucial. 

A  straightforward  treatment  of  the  analysis  of  variance  may  be 
found  in  most  texts  on  mathematical  statistics  (e.g.,  Hoel,  1947;  Rao, 
1965).  Adopting  the  notation 


JL 

ab 


I 

i.J 


ij 


Y 


•J 


i* 


t 

J-l 


U 


where  e  *»  number  of  rows  end  b  •  nueber  of  columns,  let 


(2) 


(3) 


(4) 


A  • 


-  Y,.  ♦  t  >‘ 


(5) 


j-i 


(Y 


-  Y  )' 

e  e  9 


(6) 
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b  t  Crt.  -  *..)2 

1-1  1 


(7) 


We  may  now  sake  two  hypotheses  about  the  outcome  of  the  experiment 
and,  from  Eqs.  (5)  through  (7),  form  statistics  with  which  to  test  the 
hypotheses . 


Hypothesis  1:  Removing  the  Ice  had  no  effect  the  row  means  ( Y . _ ) 
are  samples  from  the  same  distribution),  let 

« • 

Then  f  has  an  F-diatribution  with  {a  -  1,  (a  -  l)(b  -  1))  degrees  of 
freedom.  A  test  of  Hypothesis  l  is  based  on  the  following  reasoning. 

Let  Hq  represent  the  hypothesis  under  consideration  and  define 
the  conditional  probability: 

a  •  Pf{ reject  Hq1Hq  is  true) 


Next  find  the  value  F*  for  which 


F(OdC  f>t(  accept  HQ  i  «0 


irt 


true ' 


Now  if  f  %  F*  we  accept  the  hypothesis  H  (  and  if  f  *  V*  we  reject 
pe  is  found  in  any  table  of  the  F-dlstrlbuUon  as  a  fum*»  Ion  of  i  and 
the  degrees  of  freedom  (vj,v,l. 


Hypo  chest  a  II :  Perturbing  the  i  nU  la  S  condi t  f  with  a  *di  Jive  rand  a 
noise  had  no  effect  (  the  column  meat:*  tv#  •  -  r  •<  -  •  . 

distribution).  Then 


and  f  has  Jb  -  1,  (a  -  l)(b  -  l)J  degrees  of  freedom. 
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TEMPERATURE „  C.E0P0TKNTIA1,  HEIGHT.  AND  ZONAL  WIND 

The  analysis-of-variance  test  described  above  was  applied  first 
to  the  temperature,  the  geopotential  height,  and  the  zonal  (E-W)  wind 
velocity.  Our  goal  was  to  ascertain  whether  there  is  a  significant 
difference  in  these  variables  as  a  function  of  altitude  and  latitude 
when  the  arctic  sea  ice  was  removed. 

Altitude  is  here  taken  to  be  the  1000  ab,  800  mb,  and  400  mb  pres¬ 
sure  levels.  Latitude  is  usually  defined  as  an  annular  region  over 
which  the  variable  has  been  averaged.  For  instance,  we  tested  the 
value  of  the  average  temperature  between  54*N  and  70*N  latitude  at  the 
800  mb  level.  The  raw  data  going  into  the  average  were  the  values  of 
this  variable  at  all  the  grid  points  provided  by  the  simulation  lying 
within  the  specified  annular  region. 

Gradients  of  variables  have  been  specified  at  exact  latitudes  and 
have  no  meridional  averaging. 

Tabic  3  gives  all  the  pertinent  data:  the  variable,  meridional 
region,  altitude  level,  value  of  th*  variable  for  all  b  experiments 
<3  ice-in  and  3  ice-out) ,  difference  of  ice-ln  and  ice-out  mean  values, 
and  finally,  the  value  of  the  F-atatlstic  and  Its  significance  level. 

We  choose  to  accept  the  hypothesis  that  there  is  no  difference  between 
the  ice-ln  and  ice-out  results  if  the  significance  is  no  better  than 
0.05.  Therefore  we  reject  the  hypothesis  of  no  differences  if  and  only 
if  the  significance  level  y  0,05. 

WKRlPlOSAt  TRANSPORT 

In  addition  to  the  basic  variables  of  temperature,  geopotential 
height,  and  zonal  wind,  we  selected  the  stfridienat  transport  nf  sensible 
heat  and  tswscntue  for  testing.  Again,  these  derived  quantities  were  the 
30-day  mean,  serially  avenged  values.  The i r  derivation,  utilizing  the 
available  Hihta-Arakaua  data,  is  given  below: 

!iiul_ .  traps .  the  midi«Nial  sensible- heat  flux  is  given  by  I  *  ,  t  vl  , 

-1  -1 

»  *  density  el  air.  r  *  soecUie  beat  of  air  *  <>.>24  <ral  pm  deg  , 


wlie  re 


DATA  AND  :.TATliTICS  FOR  TEMPERATURE,  GE0PUTENT1AL  HEICHT,  ZONAL  UIND 
POLAR  HEAT  TRANSPORT,  AND  POLAR  MOMENTUM  TRANSPORT 


r 
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v  =  northward  wind  speed,  and  T  =  temperature.  Integrating  over  a 
rectangular  region  in  the  x  -  z  plane  gives  the  total  heat  flux,  F,  as 

/*xl  p  max 

F  a  J  I  f  dz  dx  (9) 

X0  Surface 

Upon  use  of  tne  hydrostatic  approximation  dp  *  -og  dz,  the  definition 
of  a  a-surface  dp  =  (p{  -  p^)  do,  and  the  transformation  to  polar  coor¬ 
dinates  dx  *  R  cos  dJ  ,  this  becomes 

<PS  -  Pq) vT  do  di  (10) 


e  R  cos 

f._£ - 


g 


-// 
JQ  JQ 


where  g  •  acceleration  due  to  gravity  «•  981  cm  aee~“,  p  •  surface  pres- 

s 

sure,  p0  *  pressure  at  top  of  model  atcr sphere  »  cOO  mb,  a  «•  sigma 
surface  (0(«(1),S«  radius  of  eart„  -  6.375  *  198  cm,  <t  -  latitude, 
and  i  »  longitude. 

to  utilise  the  Mitits-Arakaoa  model's  data,  which  are  given  only  at 
o  *  0.2S  and  0.75,  we  assume  that  the  mean  values  of  v  and  T  in  the  in¬ 
terval  0  S  e  £  0.5  are  given  by  the  values  at  3  *  0.25.  Similarly  the 
«Bfn  values  of  v  and  T  in  the  region  0.5  £  a  i$  1.0  arc  given  by  those 
at  e  «  0.75. 

the  vertical  integration  is  then  approximated  by 


c  R  cos  j'  y.2  * 

F  m 


2g 


i:  -  v 


)<V.25T.2S  *  V.75t.75)  d*  (U) 


With  *H  -  i-/72,  whore  72  is  the  number  of  zonal  grid  increments  in  the 
model,  the  SUB  approximating  the  meridional  heat  transport  is 


c  R  cos 
F  *  — JL 


?2g 


t  i 

Z(<".  ■  "0>1<V.2S T.»  *  V.J5T.n>1) 


(12) 


where  (...)  denotes  the  30-day  averaging  operation. 

In  practice  we  split  the  integration  shown  in  Ec>.  (12)  into  two 
part3,  the  lower  level  given  by  the  V  ^  term  and  the  upper  level 
given  by  the  V  25T  25  term*  They  may  be  added  if  the  tt  cal  heat  trans¬ 
port  is  desired.  The  reason  for  the  partitioning  is  given'  in  the  next 
section,  but  the  values  for  polar  heat  transport  for  each  level  at  70 °N 
for  the  six  cases  defined  previously  are  given  in  Table  3. 

The  F-statistics  oroduced  by  an  analysis-of-variance  test  indicated 
a  significant  difference  between  the  ice-in  and  the  ice-out  heat  trans¬ 
port  in  the  lower  level.  (We  reject  the  hypothesis  of  equality  at  the 
0.05  level) . 

There  was  no  significant  difference  in  meridional  heat  transport 
between  ice-in  and  ice-out  at  the  upper  level. 

Meridional  momentum  transport.  In  a  similar  manner  we  proceeded  to 
evaluate  the  meridional  momentum  transport  and  test  for  a  significant 
change  caused  by  the  removal  of  polar  ice. 

From  Hess  (1959) ,  the  total  poleward  momentum  transport  at  a  fixed 
latitude,  4»,  was  given  by 

x,  7. 

/I  />  max 

J  pmv  dz  dx  (13) 

*0  Surface 

where  m  »  momentum  and  v  ■  northward  wind  speed.  Making  the  substitu¬ 
tion#'  dx  ■  R  cos  <f>  dX,  m  a  (v  +  HR  cos  $)R  cos  <J>,  dp  »  -pg  dz,  we  found 

-2  2  r2v 

M  ■  - - — *  /  /  (U  +  flR  cos  <p)v  dp  dX  (14) 

*'0  *4 

r8 

where  v  ■  eastward  component  of  wind,  and  U  *»  angular  velocity  of  ro¬ 
tation  of  the  earth  •»  7.272  *  10  ^  sec 
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We  next  transform  to  o-coordinates ,  obtaining 


(p  ~  pA)  (U  +  OR  cos  4>)v  do  dX 
s  0 


(15) 


Again,  taking  the  values  at  a  =  0.25  and  a  =  0.75  as  the  mean  values  in 
the  layers  0  £  a  <  0.5  and  0.5  £  a  £  1  we  found, 


2  2  2tt 

lR  '2gS~^  Jq  (Ps  "  P0}  [(U.25  +  QR  COS  ^)v.25 
+  (U  __  +  OR  cos  (j>)v  ^  J  dX  (16) 


or 


2  2 

m  »  r*JL ..  c£3.  -1 

M  72g 


72 


i-1 


+  OR  cos  4>  ( v  2^  +  v 


]C((ps  "  P°)1tU.25v. 


,5>]) 


25  +U,75V.75 


(17) 


Table  3  shows  the  polar  momentum  transport  at  70°N  latitude.  The  ap¬ 
propriate  F-test  led  us  to  conclude  that  there  is  no  significant 
difference  in  momentum  transport  upon  removal  of  the  arctic  ice. 

Briefly  summing  up,  we  tested  a  selection  of  meteorological  vari¬ 
ables  to  ascertain  if  their  value  had  changed  significantly  after  60  days 
as  a  result  of  removing  the  arctic  se’  ice.  The  "Significance"  column 
of  Table  3  indicates  that  some  variables  changed  and  some  did  not.  Note 
that  we  have  not  estimated  a  confidence  level  for  the  magnitude  of  the 
change,  nor  have  we  discussed  steady-state  asymptotic  differences. 

Section  VIII  discusses  the  physical  implications  of  these  findings. 
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VII.  PHYSICAL  INTERPRETATION 


The  arctic  pack  ice  has  long  been  thought  to  be  a  controlling  fac¬ 
tor  in  the  climate  of  the  globe.  Budyko  (1962)  theorized  that  it  would 
require  only  small  changes  in  the  albedo  of  the  ice  to  cause  its  de¬ 
struction  within  a  few  years.  Fletcher  (1965)  pointed  out  that  the 
ice  acts  in  two  ways  to  keep  the  arctic  atmosphere  cold:  first,  it  re¬ 
flects  much  of  the  incident  solar  radiation  in  the  summer;  and,  second, 
it  prevents  heat  from  the  ocean  from  being  transferred  to  the  atmosphere 
during  the  winter.  Fletcher  further  pointed  out  that  removal  would  lead 
to  much  wanner  temperatures  over  the  Arctic  Basin  in  the  winter  but  only 
slightly  wanner  temperatures  during  the  summer.  Moreover,  an  open  Arc¬ 
tic  Ocean  would  be  a  source  of  moisture  for  the  atmosphere  that  would 
tend  to  increase  the  precipitation  in  the  Arctic. 

In  regard  to  changes  in  the  circulation  of  the  atmosphere,  Fletcher 
argued  that  a  winter  decrease  in  the  pole-to-equator  temperature  gradi¬ 
ent  would  cause  a  weakening  of  the  westerly  circulation  in  winter. 
Pogosian  (1970)  attempted  to  quantify  this  argument  by  simply  assuming 
warmer  surface  temperatures  and  constructing  maps  of  the  baric  topo¬ 
graphy  at  high  altitudes  from  simple  hydrostatic  considerations.  The 
experiment  reported  here,  carefully  designed  to  eliminate  any  differ¬ 
ences  not  attributable  to  the  removal  of  the  ice,  indicated  that  the 
model  behaves  much  in  the  fashion  that  has  been  postulated  for  the  real 
atmosphere  in  winter. 

The  most  obvious  and  luost  significant  difference  was  in  the  1000  mb 
temperature  over  the  Arctic  Ocean.  This  increase  when  the  ice  was  re¬ 
moved  is  also  significant  at  the  800  mb  level,  but  there  was  an  insig¬ 
nificant  cooling  at  the  400  mb  level.  Taken  together  they  indicated  a 
significant  steepening  of  the  lapse  rate  over  the  Arctic.  Figure  5 
shows  the  model  temperatures  on  a  tephigram  together  with  the  mean 
vertical  temperature  sounding  for  February  as  reported  by  Fletcher 
(1966b)  for  the  North  Pole  4  ice  station  during  1955-56.  The  differ¬ 
ence  in  lapse  rate  is  to  be  expected  because  with  an  open  ocean  there 
is  a  strong  heat  source  at  the  surface. 


mwmrm 
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Temperature  ( °C ) 


Fig.  5  —  Model  temperatures  on  a  tephigram  together  with  mean  vertical 
temperature  sounding  for  February  for  the  North  Pole  4  ice 
station  during  1955-56  (Fletcher,  1966b). 
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The  picture  of  the  meridional  and  vertical  distribution  of  tem¬ 
perature  difference  is  shown  in  Fig.  6.  This  figure  was  constructed 
by  forming  the  mean  monthly  zonally  averaged  temperature  for  each  lati¬ 
tude  in  the  grid  for  each  level  for  all  six  runs.  The  three  ice-in 
runs  were  then  averaged  and  subtracted  from  the  average  of  the  three 
ice-out  runs.  The  striking  feature  of  this  chart  is  the  southward  ex¬ 
tent  of  the  warning  into  sub-Arctic  latitudes.  The  vanning  was  found 
to  be  significant  in  the  latitude  band  from  54  to  70°N  at  both  the 
1000  and  800  mb  surfaces.  This  change  was  brought  about  because  the 
outbreaks  of  polar  air  which  normally  cool  the  sub-Arctic  regions  were 
not  as  cold  with  an  ice-free  Arctic  as  they  would  have  been  with  an 
ice-covered  Arctic.  This  explanation  was  further  supported  by  the 
significantly  lower  heat  transport  under  ?  -out  conditions,  across 
70*N  in  the  lower  layer  of  the  model  atmosphere.  The  changes  in  tem¬ 
perature  at  the  400  mb  level  and  south  oi  54*N  were  too  small  to  be 
ascribed  to  the  removal  of  the  ice  i  this  experiment. 

Figure  7  shows  the  changes  in  the  geopotential  height  for  the 
three  levels;  it  was  constructed  in  the  same  manner  as  Fig.  6.  The 
■oat  striking  feature  of  this  chart  is  the  lowering  of  the  geopotential 
height  in  the  lower  layers  of  the  Arctic.  If  there  were  a  strictly 
hydrostatic  change,  these  heights  would  have  been  expected  to  rise  at 
all  levels  with  very  little  rise  at  the  bottom  and  increasing  lifting 
as  altitude  increased.  The  linking  of  the  geopotential  heights  at  the 
low  altitudes  nust  have  resulted  from  a  dynamic  change  that  removed 
aase  from  the  Arctic  Region,  This  rather  small  change  In  mass  over 
the  Arctic  was  not  noticed  as  sn  Increase  in  surface  pressure  elsewhere, 
because  it  was  spread  over  e  much  larger  area  as  it  was  displaced  south¬ 
ward,  but  it  had  a  significant  effect  on  the  dynamics  of  the  circulation. 
The  two  lower  levels  were  significantly  lower  over  the  Arctic,  and  the 
two  higher  levels  were  significantly  higher  over  the  sub-Arctic.  Thin 
dynamic  redistribution  of  mass  should  have  the  effect  of  increasing  the 
weat-to-eaat  circulation  at  low  level*  in  the  Arctic,  while  decreasing 
the  west-to-east  circulation  at  high  levels  in  the  sub-Arctic.  A  test 
of  the  difference  in  geopotential  between  46  and  66*N  showed  a  signifi¬ 
cant  decrease  in  the  gradient  at  800  and  400  mb. 


ICE  EXTENT 

NORTH  LATITUDE 


Kig.  6  --  Temperature  difference* 
(*C) ;  (ice  out)  -  (ice  In*. 
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(MS) 


ICE  EXTENT 

NORTH  LATITUDE 


Fig.  ?  --  neopotentlei-lir ight  difference* 
<«) ;  (ice  out)  -  (ice  in). 
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Figure  8  presents  a  cross  section  of  the  zonal  wind  velocity  con¬ 
structed  in  the  manner  of  the  two  previous  figures.  The  decrease  of 
4  m/sec  at  the  400  mb  level  at  54°N  is  significant  and  represents  a 
confirmation  of  Fletcher's  hypothesis  for  a  weaker  circulation  in  the 
winter.  The  significant  2  m/sec  increase  at  78*N  had  not  been  pre¬ 
dicted  in  any  of  the  qualitative  studies  of  the  problem.  Although  the 
statistical  tests  failed  to  show  that  the  differences  in  the  wind  at  the 
800  mb  level  were  significant,  the  significant  changes  in  the  pressure 
gradient  responsible  for  the  increase  at  800  mb  lent  credence  to  the 
reality  of  the  change.  Another  bit  of  evidence  of  the  weakening  of 
the  sub-Arctic  upper-level  circulation  was  the  change  in  the  transport 
of  angular  momentum  into  the  Arctic.  Although  the  statistical  test 
failed  to  reach  the  5  percent  significance  level,  the  momentum  flux 
waa  cut  in  half,  indicating  a  decreaae  in  the  momentum  available  to 
be  transported  northward. 

LEVEL 

<M0) 


Fig.  8  —  Last/*t> st  wind  differences 
(a/sec);  (ice  out)  -  (ice  in). 
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VIII.  SUMMATION 


The  changes  in  the  model  behavior  between  ice-in  and  ice-out, 
which  have  been  shown  by  the  statistical  tests  to  be  the  results  of 
the  change  in  ice  conditions,  can  be  suamarized  as  follows: 

(1)  The  vertical  stability  over  the  Arctic  Basin  decreased 
markedly.  An  expected  result  of  this  decreased  stability 
would  be  an  increase  in  precipitation;  analysis  of  the 
precipitation  will  be  deferred. 

(2)  The  excess  heat  in  the  Arctic  caused  a  slight  efflux 

of  mass  from  the  area,  which  caused  lower  pressure  over 
the  Arctic  Basin. 

(3)  The  oceanic  heat  source  present  upon  removal  of  the  ice 
weakened  the  cold  core  vortex  over  the  Arctic  Basin. 

(4)  The  geopotential  heights  at  high  sltitudes  in  the  sub- 
Arctic  were  increased  by  this  change  in  the  distribution 
of  mass. 

(5)  The  mid-latitude  westerlies  were  weakened  by  the  re¬ 
distribution  of  mass. 

This  paper  has  thus  indicated  that  changes  of  the  boundary  condi¬ 
tions  in  the  model  produce  detectable  changes  in  the  output  despite  the 
breakdown  of  predictability .  We  believe  that  we  have  provided  an  ade¬ 
quate  demonstration  of  the  reaponse  of  a  primitive  aquation  model  to  a 
rather  major  change  in  the  boundary  conditions.  The  relation  between 
model  response  and  the  poaalhle  response  of  the  real  atmosphere  is  sug¬ 
gestive,  but  further  work  needs  to  be  done  to  relate  the  model  to  the 
reactions  of  the  atmosphere.  In  perusing  the  output  of  the  six  runs 
made  to  collets  this  experiment,  we  became  acutely  aware  of  the  need 
for  strict  statistical  controls.  Tha  natural  variability  of  ths  modal, 
and  certainly  tha  natural  variability  of  the  real  atmosphere,  must  al¬ 
so  be  taken  into  account  in  trying  to  assess  the  effect  of  eny  pur¬ 
poseful  or  inadvertent  changes  in  the  factors  which  are  suspected  of 
controlling  the  climate. 
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